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ABSTRACT 

We use a series of N-body/smoothed particle hydrodynamics simulations and an- 
alytic arguments to show that the presence of an effective temperature floor in the 
interstellar medium at Tp 10* K naturally explains the tendency for low-mass 
galaxies to be more spheroidal, more gas rich, and less efficient in converting baryons 
into stars than larger galaxies. The trend arises because gas pressure support becomes 
important compared to angular momentum support in small dark matter haloes. We 
suggest that dwarf galaxies with rotational velocities ~ 40 km s~^ do not originate as 
thin discs, but rather are born as thick, puffy systems. If accreted on to larger haloes, 
tenuous dwarfs of this kind will be more susceptible to gas loss or tidal transformation 
than scaled-down versions of larger spirals. For a constant temperature floor, pres- 
sure support becomes less important in large haloes, and this produces a tendency for 
massive isolated galaxies to have thinner discs and more efficient star formation than 
their less massive counterparts, as observed. 
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1 INTRODUCTION 

It is well established that small galaxies have longer star 
formation time scales than their larger cousins (e.g. Hunter 
& Gallagher 1985; van Zee 2001; Kauffmann et al. 2003; 
Brinchmann et al. 2004; Geha et al. 2006). Dwarf galaxies 
in the field have higher specific gas fractions (van Zee 2001; 
Geha et al. 2006) and are morphologically thicker (Dalcan- 
ton et al. 2004; Yoachim & Dalcanton 2006) than luminous 
late-type galaxies. It is common to discuss supernova feed- 
back or radiative feedback as a means to explain these trends 
(e.g. Dekel & Silk 1986; White & Frenk 1991; Kauffmann et 
al. 1993; Bullock et al. 2000; Mayer et al. 2001a; Somerville 
2002; Dekel & Woo 2003; Kravtsov et al. 2004; Read et al. 
2006; Stinson et al. 2006), yet the origin of the observed 
relations between galaxy properties and their total mass re- 
mains open for debate. 

In this paper we investigate the simple systematic effect 
that a gas temperature fioor, Tp , has on the star formation 
efficiencies, gas fractions, and morphologies of galaxies as 
a function of dark matter halo mass. We show that even 
a moderate effective temperature floor of Tp — 10* K, as 
might arise naturally in the presence of a photo-ionizing 
background, produces many of the general trends observed. 
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Specifically, pressure support becomes dynamically compa- 
rable to angular momentum support in small dark matter 
haloes, and this causes dwarf galaxies to be systematically 
puffier and less efficient in converting gas into stars than 
their larger counterparts. 

In what follows we explore a range of values Tf = 
(1.5 — 5) X 10* K, in order to qualitatively mimic real phys- 
ical effects that may act to pressurise gas within galaxies. 
These may include inefficient cooling, heating by an internal 
or external ultraviolet (UV) background, supernova feed- 
back, turbulent pressure, and cosmic-ray heating, among 
others. Of course, within a fully realistic galactic interstellar 
medium (ISM) we expect molecular cooling to produce cold 
clouds embedded within a warm, pressurised medium (e.g. 
Mckee & Ostriker 1977; Yepes et al. 1997; Springel & Hern- 
quist 2003; Robertson et al. 2004, 2006) Internal and exter- 
nal heating sources and cooling within the medium should 
produce a quasi-stable system with an effective temperature 
that acts to stabilise the galaxy. Unfortunately, the nature of 
the coupling between various energy sources and the back- 
ground ISM is poorly understood. Rather than attempt to 
model these processes, we use Tf as a phenomenological 
proxy. 

Our approach is perhaps most directly relevant to the 
situation of an ionizing background field that acts to prevent 
cooling below ~ 10* K in warm galactic gas. It is well known 
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that a UV background can suppress galaxy formation alto- 
gether in very small galaxies < 30kms~^ (Efstathiou 1992; 
Thoul & Weinberg 1996; Quinn et al. 1996; Gnedin 2000; 
Hoeft et al. 2006; Grain et al. 2007) and that this effect can 
have important implications for the baryonic mass fractions 
and overall abundance of the smallest galaxies in the uni- 
verse (e.g. Bullock et al. 2000; Benson et al. 2002; Strigari 
et al. 2007). In a recent examination, Hoeft et al. (2006) 
used cosmological simulations with cooling and star forma- 
tion to show that gas at typical galactic densities will have 
equilibrium temperatures ~ (1 — 3) x lO** K, and that this 
temperature is fairly insensitive to the normalisation of the 
UV flux. This heating strongly suppresses the baryon frac- 
tion in haloes smaller than ~ 20 kms~^ in their simulations. 
Our investigations focus on galaxies that are just above this 
scale. Specifically, we explore the morphology and star for- 
mation efficiency of galaxies that are large enough to accrete 
warm gas but small enough to be dynamically affected by 
warm gas pressure. 

We note that the effect of a finite temperature floor 
on dwarf galaxy formation was discussed in a semi-analytic 
context by Kravtsov et al. (2004), who used the idea to mo- 
tivate models for low star formation rates in small galaxies; 
and by Tassis et al. (2006), who used a similar model to in- 
vestigate the gas fraction and mass-metallicity relationships 
in dwarf galaxies. Taylor & Webster (2005) studied star for- 
mation within equilibrium dwarf galaxies by modeling H2 
cooling within a thermally-balanced ~ 10* K medium and 
derived lower limits on self-regulated star formation rates in 
dwarfs in this context. 

In the next section we present an analytic investigation 
into the importance of baryonic pressure support compared 
to angular momentum support as a function of virial mass 
and gas temperature and show that pressure support should 
become dynamically important in dwarf-galaxy haloes. In §3 
we use the N-body/smoothed particle hydrodynamics (SPH) 
code Gasoline (Wadsley et al. 2004) to explore galaxy for- 
mation with a range of temperature floors and halo masses. 
We present results on galaxy disc thickness, gas fractions, 
and star formation rates as a function of galaxy circular ve- 
locity. We reserve §4 for discussion and §5 for conclusions. 



2 ANALYTICAL EXPECTATIONS 

The standard analytic approach for calculating galaxy sizes 
and morphologies within dark matter haloes assumes that 
the gas cools to a temperature well below the halo virial 
temperature, Tg ^ Tv. As a result, the thermal pressure 
support in the gas is small compared to its angular mo- 
mentum support. A thin disc of star- forming material is 
the natural outcome (Fall & Efstathiou 1980, Blumenthal 
et al. 1986). This thin-disc configuration is taken to be the 
starting point for galaxy formation (and star formation) in 
most models, including those that model dwarf galaxies (e.g. 
Kauffmann et al. 1993; Somerville & Primack 1999; Benson 
et al. 2002; Somerville 2002; Mastropietro et al. 2005; Mayer 
et al. 2006; Gnedin et al. 2006; Button et al. 2007). Here we 
stress that the thin disc approximation should break down 
in small haloes where TV Tg ~ 10* K. In these cases, the 
pressure support radius becomes comparable to the angular 
momentum support radius and we expect a thick morphol- 



ogy. What follows is a simple analytic investigation aimed 
at quantifying the halo mass scale of relevance. A more rig- 
orous numerical approach is given in the next section. 

Consider dark matter haloes of mass Mv with virial 
radu definecQ such that ~ 113 kpc (Mv/lO^Mo)!/^ We 
assume that halo density profiles are well approximated by 
the NFW fit (Navarro, Frenk, & White 1996): 



where x = r/vs and the scale radius, Vs, is determined 
from the halo concentration parameter Cv = Rv/vs- We 
adopt the relation Cv = 10 (Afv/lO" Mo)"""**^, which is 
appropriate for a ag = 0.75 AGDM cosmology (Bullock 
et al. 2001a; Maccio et al. 2007). Given and Cv, the 
circular velocity curve, Vc{r) = GM{r)/r, is determined 
by the integrated mass profile. For our adopted relation, 
the circular velocity peaks at a maximum value Vmax — 
71kms-i(Mv/10" Mq)°-'\ 

In the standard thin-disc scenario, the gas obtains 
specific angular momentum that is similar to that of the 
dark matter, which is often characterised by a dimen- 
sionless spin parameter (Peebles 1969) defined as A = 
j\E\^^^ M~^''^, where G is Newton's constant and j 
and E are the speciflc angular momentum and energy of 
the halo, respectively. Simulated CDM haloes typically have 
A ~ 0.03 with a 90 % spread between 0.01 — 0.1 (e.g. Barnes 
& Efstathiou 1987; Bullock et al. 2001b; Maccio et al. 2007). 

It is straightforward to show that if the gas cools and 
contracts without angular momentum loss to form a thin, 
angular momentum supported exponential disc, the disc 
scale radius is given by (Mo, Mao, & White 1998; hereafter 
MMW) 

i?d = Ai?haio/(c, A,md). (2) 

Here we assume that the gas falls in from a radius -Rhaio, 
deflned to be either the virial radius, i?v, or the "cooling 
radius", Rc, depending on which one is smaller -Rhaio = 
min(_Rc, -Rv). By introducing the cooling radius (White & 
Frenk 1991) we account for the expectation that hot gas 
in the outskirts of massive haloes will not have had time 
to cool since the halo formed. For simplicity, we adop10 
-Rc ~ 129 kpc V^2o''*i where V120 is the halo maximum circu- 
lar velocity in units of 120 km s~^. The function / (~ 1) in 
Equation [2] contains information on the halo profile shape 
and contraction from baryonic infall, and depends on the 
initial halo concentration c = Rhaio /Ra and the disc mass, 
rud, in units of the total mass within iihaio. 

The three solid lines labelled with A values in Figure [T] 
show Rd calculated in the standard thin-disc framework as a 
function of the (initial) halo Vmax value. From top to bottom 
the lines assume spin parameters A = 0.1, 0.03, and 0.01. We 
have used the fitting formula from MMW /(c, A, m^) with 
TUd = 0.1. We plot the scale radius as a function of the initial, 
uncontracted halo Knax in order to facilitate the following 
comparison. 

^ In our definition, the virial radius contains an average mass 
density equal to 360 times the matter density of the universe and 
we adopt Qm = 1 — f^A = 0.27 and with h = 0.7. 
^ This approximation was given by Mailer & Bullock (2004) for 
haloes with V120 ^ 1- Below this scale, the virial radius sets -Rhaio- 
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Table 1. Simulated Galaxies 
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(7) 
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Name 




Cv 
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Ri 




V2.2 




n —11 
[kms J 




[10^" Mq] 


r-1 n4 T^i 
[10 KJ 


[kpc] 


n 1 
[kpcj 


n —11 
[kms J 


D24 


24 


14 


0.29 


1.5 


0.37 


0.08 


37 


D24_g 


24 


14 


0.29 


1.5 


0.35 


0.19 


35 


D28 


28 


13 


0.47 


1.5 


0.55 


0.09 


43 


D28_g 


28 


13 


0.47 


1.5 


0.56 


0.18 


41 


D41 


41 


12 


1.6 


1.5 


0.92 


0.10 


65 


D41.g 


41 


12 


1.6 


1.5 


1.02 


0.16 


65 


D41.W 


41 


12 


1.6 


3.0 


0.61 


0.13 


64 


D53 


53 


11 


3.8 


1.5 


1.32 


0.12 


84 


D53.g 


53 


11 


3.8 


1.5 


1.34 


0.15 


84 


D53.W 


53 


11 


3.8 


3.0 


1.05 


0.16 


82 


D53.H 


53 


11 


3.8 


5.0 


0.77 


0.18 


84 


G74 


74 


10 


11.0 


1.5 


1.6 


0.14t 


115 


G74.g 


74 


10 


11.0 


1.5 


2.2* 


0.03* 


92 


G148 


148 


8 


100. 


1.5 


1.88 


0.22t 


209 


G148.g 


148 


8 


100. 


1.5 


1.9* 


0.02* 


217 


G168.gHRLS 


168 


11 


115. 


3.0 


1.94 


0.13 


213 



(1) Names labelled "_g" refer to pure gas runs without star formation whereas names with "_W" and "_H" signify "warm" or "hot" 
temperature floors and "HRLS" refers to the Milky Way model described in Kaufmann et al. (2007). (2-4) The listed values for Vmax, 
Cv and Mv are initial halo parameters. (5) Tp is the imposed temperature floor. (6-8) The last three columns list parameters measured 
in the flnal galaxy, where Rd and Zd are the exponential scale length and scale height, respectively, and V2.2 is the rotational velocity in 
the gas measured at 2.2 R^. Disc scale heights marked with (f) are artificially large as a result of numerical heating. Gas disc 
parameters marked with (*) are derived from discs which have become Toomre unstable; thin, but spatially irregular. 

Consider now the galaxy radius that would result from 
pressure support in an idealised, spherically symmetric sys- 
tem consisting of a gravitationally sub-dominant gas. As- 
sume that this gas has no angular momentum but reaches 
a temperature of Tg = Tp within an extended NFW halo of 
virial temperature Tv. We emphasize that we will calculate 
the pressure support radius for "cold" gas at Tg = Tp, the 
lowest temperature the gas can reach within our assump- 
tions. For this approximate calculation, we neglect baryonic 
contraction. We will define the virial temperature in analogy 
with an isothermal gas of temperature T and its associated 
speed of sound Cs (e.g. Mailer & Bullock 2004): 

Specifically, Tv is set by using Knax/\/2 for Cs in Equation 

El 

The equilibrium gas profile, pgij-), will be set by a com- 
petition between the isothermal gas pressure, Pg = CgPg, 
and the gravitational potential. If we assume that the grav- 
itational force is dominated by the NFW halo potential, the 
hydrostatic force balance equation 

^dpg_ ^ ~V^{r) ^ 
pg dr r ' 

can be rewritten in terms of the dimensionless radial param- 
eter X = r/vs as 

7gt - 

Here, rj = T^/Tp parameterises the relative strength of 
the halo gravity and the thermal pressure of the gas and 
h{x) = 9.26a;"^ [ln(l -l-a;) - 2:/(l -l-a;)]. Note that in the 
limit of large 77 (Tv ^ Tp), the gas profile will be centrally 
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Figure 1. Galaxy radius R* versus halo Vmax for various models 
for galaxy assembly. _R* is defined as the radius that contains 
26% of the galaxy mass and is equivalent to the disc scale radius 
for a thin exponential disc. The solid lines that increase with 
decreasing halo size show R* for a pressure-supported isothermal 
gas with temperature floors as labelled. The straight lines that 
break at Vmax ~ 90kms~^ show R* for an angular momentum 
supported thin disc for three different halo spin parameters as 
labelled. The break occurs when the infall radius transitions from 
the virial radius (Rv) to the cooling radius (-Rcool)- Note that the 
pressure support radius becomes comparable to (and even larger 
than) the angular momentum support radius in small galaxies. 
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concentrated with a large negative derivative, and a negligi- 
ble pressure-support radius. More generally, solving Equa- 
tion [5] for pg yields 



Paix) = poexp 



-9.26 »7 {l 



ln[l + x] 



(6) 



where the normalisation parameter po sets the gas density 
at a; = 0. It is clear that in small haloes with 77 ~ 1, the gas 
profile can extend to a; ~ f or r ~ _Rv/c ~ O.f-Rv, which is 
comparable in size to the angular momentum support radius 
(Eq.0. 

The thick-solid and dotted lines in Figure [T] show a 
more explicit comparison for various gas floor temperatures: 
Tf — f, f.5, and 5 x 10'' K. We have used Equation [S] to 
compute the radius K* that encloses 26% of the pressure- 
supported galaxy mass as a function of halo Knax- This ra- 
dius is analogous to the scale radius for an exponential disc, 
which contains 26% of the rotationally supported disc mass. 
While relatively unimportant in large Milky- Way-size haloes 
(Knax ~ 200kms~^), we see that pressure support should 
dominate in shaping galaxy morphologies and gas distribu- 
tions in small haloes. The expectation is that small galaxies 
will be intrinsically puffier than large galaxies, even in the 
absence of environmental influences. 

The scale where pressure becomes important compared 
to rotation will naturally depend on the temperature floor 
of the gas and on the intrinsic spin. For Tf — 1.5 x 10* K 
and A = 0.03, we expect the effect to become very important 
for dwarf-size haloes with Vmax ^ 35kms~^. If the temper- 
ature floor is high, Tp = 5 x 10* K, then the effect could be 
important even in ~ lOOkms"^ haloes if they have inhabit 
the low-spin tail of the distribution, A = 0.01. If galaxies 
form with the range of spins expected (A ~ 0.01 — 0.1) then 
we would predict a range of morphologies (from puffy to 
thin discs) at fixed Vmax as long as the temperature floor is 
roughly the same from galaxy to galaxy. Though we do not 
explore galaxy formation in very small halos (< 20km s~^) 
in the rest of this paper, it is interesting to note that we 
expect the morphologies of the smallest objects to be essen- 
tially spheroidal, with initially extended gas profiles. The 
stellar sizes of these objects will likely be much smaller than 
the gas extent, as high-densities will be required for star 
formation. This might provide an explanation for why the 
smallest galaxies (dwarf spheroidals) are always dispersion 
supported systems. 

The simple, spherical, model we have just explored was 
primarily designed to guide expectations. Stronger results 
are presented in the next section, where we use 3D hydro- 
dynamical simulations to investigate the effect of a reason- 
able gas temperature floor on morphologies, gas fractions, 
and star formation rates in small galaxy haloes. 



3 SPH SIMULATIONS 

We use the parallel TreeSPH code Gasoline (Wadsley et al. 
2004) , which is an extension of the pure N-Body gravity code 
Pkdgrav developed by Stadel (2001). It includes artificial 
viscosity using the shear reduced version (Balsara 1995) of 
the standard Monaghan (1992) implementation. Gasoline 
uses a spline kernel with compact support for the softening of 
the gravitational and SPH quantities. The energy equation 



is solved using the asymmetric formulation, which is shown 
to yield very similar results compared to the entropy con- 
serving formulation but conserves energy better (Wadsley et 
al. 2004). The code includes radiative cooling for a primor- 
dial mixture of helium and (atomic) hydrogen. Because of 
the lack of molecular cooling and metals, the efficiency of 
our cooling functions drops rapidly below 10* K. The lack 
of molecular cooling is unimportant in our investigation be- 
cause we enforce temperature floors Tf 1.5 x 10* K. 

We investigate runs with and without star formation. 
The adopted star formation recipe is similar to that de- 
scribed in Katz (1992). Specifically, a gas particle may spawn 
star particles if i) it is in an over-dense region; ii) it is cool, 
with T — Tf; and iii) it has a density greater than a critical 
threshold, pg > psF- In practice, the critical star formation 
density is the most important parameter. In our primary 
simulations we use psF ~ 2.5 x lO^Mokpc""^, but explore 
a case with psF increased by a factor of 100 in §4. Once a 
gas particle is eligible for spawning stars, it does so based 
on a probability distribution function with a star forma- 
tion efficiency factor c* that is tuned to match the Kenni- 
cutt (1998) Schmidt Law for the Milky Way and M33-size 
disc described in Kaufmann et al. (2007). The mass of gas 
particles decreases gradually as they spawn more star par- 
ticles. After its mass has decreased below 10% of its initial 
value the gas particle is removed and its mass is re-allocated 
among the neighboring gas particles. Up to six star particles 
are then created for each gas particle in the disc. We note 
that Stinson et al. (2006) have implemented a similar star 
formation recipe, although they include an allowance for su- 
pernova feedback effects using a subgrid, multi-phase model 
based on blast waves. 



3.1 Initial conditions 

We simulate 15 isolated systems with masses spanning the 
scale of dwarf galaxies to large spirals, with initial maximum 
circular velocities that range from V^max = 24 to 148kms~^. 
We also present results from an older simulation of a Milky- 
Way size galaxy (Knax = 168kms~^) that was originally 
discussed in Kaufmann et al. (2007). Haloes are initialised 
as spherical equilibrium NFW prof iles using the methods 
outlined in Kazantzidis et al.l (|2004f ) and we use the mass- 
concentration relationship discussed in §2 to set the profile 
parameters. Table 1 lists the specific parameters used in each 
simulation and provides a reference name for each run. 

We initialise a fraction of the total halo mass, /b — 0.1, 
as a hot baryonic component with the same radial distribu- 
tion as the dark matter and impose a temperature profile 
such that the gas is initially in hydrostatic equilibrium with 
an adiabatic equation of state. For all of our fiducial mod- 
els we choose Xg = 0.03 for our gas spin parameter, defined 
in analogy with the halo spin as Xg = jg\E\^^'^G~^My^^'^. 
Here, jg is the average specific angular momentum of the 
gas, E and Afv are the total energy and mass of the halo. 

The specific angular momentum distribution of the gas 
is assumed to scale linearly with the cylindrical distance 
from the angular-momentum axis of the halo, j oc r^". 
This choice is consistent with values found for dark mat- 
ter haloes within cosmological N-body simulations (Bullock 
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et al. 2001b). For simplicity, we initialise the dark matter 
particles with no net angular momentum. 

The hot gaseous halo is sampled with Ng = 10^ particles 
and the dark matter halo with TVdm = 2 x 10^ particles. The 
force resolution is set to be a fixed fraction of the simulated 
halo virial radius, /res = 0.002_Rv. These choices correspond 
to cases where numerical losses of angular momentum be- 
come small (Kaufmann et al. 2007). A detailed description 
of our initialisation method and results on the evolution of 
a Milky Way-size galaxy is presented in Kaufmann et al. 
(2007). 

As seen in Table 1, our new simulations sample six 
halo circular velocity scales spanning those of dwarfs to 
large spirals. For each Vmax, we ran a case with and with- 
out star formation and used a conservatively low temper- 
ature floor Tf = 1.5 X 10'' K. We performed two addi- 
tional Vmax ~ 53kms~^ simulations with a "warm" and 
"hot" temperature floors at Tf ~ 3 and 5 x lO' K, respec- 
tively, and simulated a second V^ax = 41kms~^ case with 
Tf = 3 x 10* K. Higher temperature floors were not explored 
in the smaller haloes because galaxy formation is suppressed 
all together if Tf ~ Tvlfl 

The galaxies were evolved for 5 Gyr, but we flnd that the 
global results stabilise after 3 Gyr (see Figure 6 below) . The 
two largest "gas-only" runs without star formation (G74_g 
and G148_g) produced gas discs that were unstable to their 
own self-gravity and became clumpy. This is not too sur- 
prising given the low Tf adopted, and is consistent with 
previous claims that substantial heating of the interstellar 
medium is needed to stabilise disc galaxies (e.g. Robertson 
et al. 2004) . While we have listed radially and vertically av- 
erage "disc" properties for these unstable cases in Table 1, 
we do not include these systems in relevant figures below. 
In order to present results for a large, stable, pure gas disc, 
we use the Vmax ~ 168 kms~^ simulation from Kaufmann et 
al. (2007) with Tf = 3 x lO" K. This system had a slightly 
larger spin parameter than the rest of our runs (Ag = 0.038). 

We note that our final galaxies have larger maximum 
circular velocity scales than their initial haloes because of 
the effects of baryonic contraction (compare the second and 
last columns in Table 1). For example, our D41 series pro- 
duces galaxies that are comparable in rotation speed to 
dwarf irregular galaxies like the Large Magellanic Clouds 
at ~ 60kms~'^. Our smallest (D24) runs produce galaxies 
that are large enough (~ 35kms~^) to be included in the 
Geha et al. (2006) sample of SDSS dwarfs. Our larger galax- 
ies, G74 and G148, produce systems that are comparable to 
M33 (~ 100 km s-') and the Milky Way (~ 200 km s"'). 

Finally, we address our initial conditions in light of the 
idea of "cold flows" put forward by Keres et al. (2005) and 
Birnboim & Dekel (2003). These authors find that gas is 
not shock-heated to the virial temperature in small haloes 
Mv < IO'^^Mq, but rather is accreted as "cold" material, 
with Tg < 10^ K. Our models are not strongly at odds with 
this picture. Specifically, the gas within our small haloes 



^ Of course, in a more complicated scenario, feedback from star 
formation may be the source of the a high effective gas tempera- 
ture, in which case a dwarf galaxy could form, and subsequently 
lose its gas once it is heated to Tg ~ Tv (see Stinson et al. 2007). 




Figure 2. Edge-on views of a simulated dwarf galaxy (D24, up- 
per) and a more massive galaxy (D53, lower). The gray-scale 
maps the projected stellar density. The disc of the larger galaxy 
is clearly thinner than the disc of the dwarf. The vertical bars 
indicate twice the softening length used in the simulations. 
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Figure 3. Galaxy disc "thinness" as a function of circular ve- 
locity for our star formation runs. Here V2.2 is the rotational 
velocity of the gas measured at 2.2_R;j for each disc. Symbol types 
correspond to different temperature floors as indicated. Discs are 
thicker in smaller haloes and for larger temperature floors. The 
arrows indicate that the stars in the largest galaxies have been 
artificially thickened by numerical heating. 



cools very quickly to the temperature floor and indeed falls 
into the central region in its "cold" phase. 



3.2 Results 1: morphological trends 

The upper and lower panels of Figure [2] illustrate the flnal 
projected stellar density for the D24 and D53 runs. It is 
evident that even the mild temperature floor, Tp = 1.5 x 10'* 
K has resulted in a very thick disc for the small, ~ 35 km s~ ^ 
galaxy, while the larger system, ~ 85kms~'^, is closer to a 
standard thin disc. 
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Figure 4. The disc "thinness" as a function of galaxy rotational 
velocity at 2.2 Rd for our pure gas runs. Symbol types correspond 
to different temperature floors as indicated. The morphological 
trend is quite similar to that obtained for our runs with star 
formation (Figure [Sjl. 

We note that both of the systems in Figure |2] have 
formed substantial central bulge components. These central 
stellar mass concentrations are likely an artifact of our sim- 
ple initial conditions, which assume a centrally-concentrated 
NFW profile for the hot gas. We return to this potential 
shortcoming in §4. 

In order to provide a more quantitative comparison be- 
tween runs, we have estimated a disc scale length, Rd, and 
scale height Zd, for each galaxy. We have explored several 
methods for quantifying Rd and Zd and find that our overall 
results change very little between methods. Because of the 
large bulge component, our galaxies are not well described 
by a single exponential surface density profile. Instead of 
adopting a degenerate two-component fit, we define Rd by 
simply measuring the radius where the face-on stellar (or 
gas) surface density drops to 0.03 (or e"'' *^) of its central 
value, and define this radius to be equal to 3.45-Rd- This 
method is generally insensitive to the details of the func- 
tional form of the final disc and probes a radial range that 
is well sampled with star (or gas) particles. The vertical scale 
height, Zd is determined by fitting an exponential profile to 
the projected, edge-on surface density profile at a projected 
radius equal to Rd (this avoids the bulge region). The fit 
is typically good out to vertical scales as large as z ~ 3 
above the disc. The measured values for each simulation are 
listed in Table 1. 

Figure |3] presents the galaxy "thinness" ratio {Rd/zd) 
as a function of the (edge-on) gas rotational velocity, V2.2, 
measured at 2.2_Rd. The points indicated by stars correspond 
to the Tf = 1.5 X 10* K cases, while the solid hexagons and 
open diamond show Tf — 3 and 5 x 10* K, respectively. It is 
clear that at a fixed Tf small galaxies are naturally thicker 
than large galaxies. This simple prescription reproduces well 
the general trend found by Yoachim & Dalcanton (2006) for 
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Figure 5. Star formation efficiency (see text) as a function of the 
galaxy rotational velocity. Symbol types correspond to different 
temperature floors as indicated. 

edge-on galaxy thickness as a function of circular velocity 
(e.g. their Figure 5). 

We note that the two largest galaxies in Figure|3]are ar- 
tificially thickened (as indicated by the arrows). Both force 
softening and two-body heating are important in these sys- 
tems. The force softening acts as an artificial pressure at 
small scales (Bate & Burkert 1997) and this effect becomes 
important in the two most massive galaxies G74 and (espe- 
cially) G148, where the vertical extent of the disc becomes 
comparable to softening parameter (/res — 0.25 and 0.5 kpc, 
respectively). Furthermore, because the massive galaxies are 
more efficient in spawning star particles, the ratio of the 
average star particle mass to dark matter particle mass is 
higher in the larger galaxies (the ratio is one order of mag- 
nitude higher in G148 compared to D24). We therefore ex- 
pect that our measured Zd values are over-estimates of the 
vertical scale heights compared to what would have been 
achieved with higher resolution simulations (which would be 
numerically very expensive). These numerical problems act 
preferentially to thicken the larger galaxies. Therefore, the 
general trend with decreasing thickness see between ~ 40 
and lOOkms"^ should be robust. 

Figure |4] shows that the same morphological trend is 
seen in our our pure gas runs without star formation. We do 
not plot the massive G74_g and G143_g galaxies. As men- 
tioned above, the gas became so cold in these runs {Tf = 
1.5 X 10* K) that the final discs fragmented into thin, clumpy, 
irregular systems. The open hexagon at V2.2 ~ 213kms~^ 
is simulation G168_gHRLS, where HRLS refers to the Milky 
Way model described in Kaufmann et al. (2007). This sys- 
tem is hot enough to be stable, with Tp = 3 x 10* K, and 
ends up as very thin disc. Overall, the agreement between 
our pure gas runs and those with star formation suggests 
that the correlation between disc thinness and circular veloc- 
ity should hold, and is largely independent of uncertainties 
associated with star formation. 
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3.3 Results 2: gas fractions and star formation 

In the previous section we sliowed tliat galaxies formed 
witiiin small lialoes tend to be tliicker tiian tiiose formed 
witiiin large lialoes. Figure [5] sliows tliat tlie star formation 
efficiency (SFE) in our simulated galaxies also varies as a 
function of ¥2.2- We define SFE = m«/mg, where rh* is the 
star formation rate and rrig is the gas associated with the 
galaxy. Specifically, rrig is defined to be the mass of gas that 
is both cold (T — Tf) and no longer infalling. It is evident 
from Figure [5] that our dwarf galaxies are less efficient in 
turning gas into stars than are larger galaxies, as expected. 
Moreover, at fixed circular velocity, the efficiency is reduced 
for higher ISM temperature floors. 

Another observationally-oriented measure of the effi- 
ciency of star formation is the cool gas fraction, mg/(mg + 
m*). The upper and lower panels of Figure [S] show the evo- 
lution of the cool gas with time in our galaxies. We find that 
the cool gas fraction approaches a constant after ~ 3 Gyr. 
This corresponds to the time when the infall of new cool 
gas reaches an equilibrium with the rate that gas is being 
converted into stars in the galaxy. The upper panel shows 
the evolution for galaxies with different (final) circular ve- 
locities (V2.2) at a fixed ISM temperature Tp = 1.5 x 10* K. 
The lower panel shows our D53 series for three values of Tf- 
Larger systems end up with lower gas fractions, as do sys- 
tems with decreasing ISM temperatures. Figure [7] shows the 
same data sliced at a fixed time (3 Gyr) plotted as a func- 
tion of V2.2- Clearly, the gas fractions are higher in smaller 
galaxies, as expected. 



4 DISCUSSION 

Geha et al. (2006) looked at a sample of 101 extremely low 
luminosity dwarf galaxies selected from the Sloan Digital 
Sky Survey, and found a trend for dwarfs to be systemati- 
cally much more gas rich than giants (e.g. Geha et a. 2006, 
Figure 3) . The results presented in the previous section (Fig- 
ureO show encouraging agreement with the observed trend. 
However, Geha et al. find an average gas fraction in dwarfs 
of ~ 0.6, and our simulated dwarfs have gas fractions that 
are lower ~ 0.3. The comparison may be even worse than it 
appears because the observations constrain only the neutral 
infraction. At T ~ 10* K, however, the difference between 
neutral and total fraction is expected to be small. While we 
regard the predicted trend between gas fraction and veloc- 
ity scale as the most robust aspect of this work, it is worth 
investigating whether simple adjustments might bring our 
results into closer agreement with the dataQ 

An obvious problem with our simulations is that our 
galaxies all form pronounced bulges (see Figure 2). As dis- 
cussed elsewhere, the bulges may be artifacts of our simplis- 
tic initial conditions, which assume pure NFW profiles for 
the hot gas haloes (Kaufmann et al., in preparation; see also 
Hansen & Sommer-Larsen 2006; Mastropietro et al. 2005). 

* Note, liowever, that these observations possibly underestimate 
the stellar masses because contributions from extended, low sur- 
face brightness stellar populations are likely to be missed (M. 
Blanton, private communication; or see Roberts & Haynes 1994; 
van Zee 2001). 



a 



0.2 







;r, 1 1 1 1 


V ~ 


1 1 1 1 1 

37 km S-' - 


' \ \ \ 

- W \ \ 

- w \ \ 


V ~ 


43 km S-' 


V ~ 


65 km s"' - 


V ~ 


84 km S-' - 


- v\ ^ 


V ~ 


115 km S-' - 


\ \ \ ■■■■■ 

Y ^. 

\ \ 


V ~ 


209 km S-' - 








. . . 1 . 




1 . . . . 



I I I I I I I I I I I I I I I L 

12 3 

time [Gyr] 




= 50000K 
= 30000K 
= 15000K 



0.2 - 





12 3 4 

time [Gyr] 

Figure 6. The fraction of cool gas in galaxies as a function of 
time. The upper panel shows runs with different halo masses (V2.2 
Hsted) and a fixed ISM temperature floor Tf = 1.5 X 10* K. The 
lower panel shows a series of models with the same halo mass (the 
D53 series) but different Tf values. 

In Figure [SI we have attempted to account for this obvious 
shortcoming by plotting the disc gas fraction as a function 
of circular velocity for the same set of galaxies shown in 
Figure [T] Here, we have removed the central stellar compo- 
nent in computing the disc gas fraction, mg/(mg -I- m.disc)- 
Specifically, we do not include central stellar material with a 
spherically-averaged density larger than 2.3 x 10^ Mq kpc~^, 
and this eliminates the bulge-like nuclear structures in our 
galaxies. As seen by the stars, hexagons, and open diamond 
in Figure 8, these "disc" gas fractions are more in line with 
observations. 

Of course, another clear uncertainty in any galaxy for- 
mation simulation is star formation. The most important 
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parameter in our prescription is the density threshold for 
star formation, psF- Observationally, we can constrain only 
the relationship between the projected gas density and star 
formation rate (Kennicutt 1998), however, because the den- 
sity threshold in our prescription is three-dimensional, we 
are left with the freedom to explore its parameter space 
The squares with crosses in Figure |8] show the result of two 
runs (D24_d and G148_d) with psF set at 100 times our fidu- 
cial value (to 2.5 x 10® Mq kpc~^). Note that in order to 
more directly compare with the fiducial runs, we have not 
excluded bulge stars in these points. As might be expected, 
the increased threshold produced a much higher gas frac- 
tion for the dwarf galaxy compared to the standard case 
(Figure 7). It also produced a smaller disc than the fiducial 
run, but a similar axis ratio. Unlike the small galaxy, the 
gas fraction in the G148_d system is quite similar to that 
in the fiducial run (Figure 7) because the gas was able to 
become quite dense and form stars. Like the fiducial Milky 
Way-size galaxy, this system also sits on the Kennicutt rela- 
tion. We note, however, that unlike the fiducial case, G148_d 
produced ~ 10 very dense star clusters in the final disc. Nu- 
merical scattering off of these clusters dramatically affected 
the star particle orbits, and artificially increased the scale 
height of the final disc. 

Elmegreen & Parravano (1994), Blitz & Rosolowski 
(2004, 2006) and Wong & Blitz (2002) pointed out the rela- 
tion between the gas pressure in the midplane of the galaxy 
and the ratio of H2 versus H I , and related also the star for- 
mation efficiency. The gas pressure in the midplane P/fcs in 
our simulations has values ranging from 4 X 10* to 8 X lO'^ 
cm~'^ K in agreement with the findings of Blitz & Rosolowski 
(2006), but does not show a strong evolution with galaxy 



^ We refer the reader to Kravtsov (2003) for an interesting theo- 
retical discussion on the origin of the Schmidt-Kennicutt relation. 
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Figure 8. Stars, hexagons, and open diamond show the cool gas 
fraction in the disc for the fiducial series of simulated galaxies 
(bulge stars have been removed from this measurement, see text 
for details). As before, the symbol types indicate the temperature 
floor assumed. The crossed squares show the galaxy gas fraction 
in runs with an increased star formation density threshold pgp — > 
lOOpsp. For those runs the bulge stars have been included in the 
measurements. 



mass. In the simulations the pressure stabilises the gas well 
above the midplane for the dwarfs, but not for the more mas- 
sive galaxies, making the midplane pressure less meaningful 
for the dwarf galaxies. 

While it is not the aim of this paper to reproduce the 
observations in detail, we conclude that there are several 
physically plausible effects that can eventually lead to a 
more complete understanding of the gas fractions in small 
systems. Overall, the agreement between the predicted and 
observed trends is quite encouraging. 



5 SUMMARY AND CONCLUSIONS 

We have used SPH simulations and an analytic discussion to 
argue that many of the observed changes in galaxy proper- 
ties as a function of their rotation speed arise naturally be- 
cause of the increased importance of ISM pressure in small 
haloes. Our main results may be summarised as follows: 

• An effective ISM temperature floor at ~ 10* K produces 
a pressure support radius that is comparable to the angular 
momentum support radius in dwarf galaxies, ~ 40kms^^. 
This suggests that most small galaxies are not formed as 
thin discs, but rather are born as thick, puffy systems. 

• For a constant temperature floor, pressure support be- 
comes less important in large haloes, and this naturally pro- 
duces a tendency for massive, isolated galaxies to have thin- 
ner discs than their less massive counterparts, as observed. 

• The morphological trend produces related trends in star 
formation efficiency; dwarf galaxies are predicted to have 
longer star formation time scales than larger galaxies. Simi- 



larly, galaxy gas fractions decrease with circular velocity, as 
observed. 

The expected morphological trend seems to be fairly 
independent of star formation details (compare Figures O 
and|4|. While relations of this kind are difficult to quantify 
observationally for a large number of galaxies, Yoachim & 
Dalcanton (2006) used a sample of 34 late-type, edge-on, 
undisturbed disc galaxies to show that more massive galaxies 
are generally thinner than less massive galaxies (see their 
Figure 5). They find radial to vertical axis-ratios for dwarf 
galaxies as low as ~ 3, in agreement with results presented 
in our Figure |3l 

More clues to the nature of galaxy formation on small 
scales can be gained from the population of dwarf spheroidal 
galaxies (dSphs) . Unlike equally-faint dwarf irregular galax- 
ies (dl), dSphs are gas poor. The fact that gas-poor dwarfs 
are exclusively found in the proximity of a luminous neigh- 
bor (Geha et al. 2006) encourages the notion that tidal forces 
and ram pressure stripping act to transform dl type galaxies 
into dSphs (Gunn & Gott 1972; Lin & Faber 1983; Moore 
& Davis 1994; Mastropietro et al. 2005; Mayer et al. 2QQla, 
2001b, 2002, 2006, 2007). Typically, models aimed at testing 
this transformation hypothesis initialise a thin disc within 
a small dark matter halo and investigate how tides or ram- 
pressure affect the galaxy as it falls into a larger host. 

Mayer et al. (2006, 2007) used Smoothed Particle Ify- 
drodynamics (SPH) simulations to show that the combined 
effects of tides and ram pressure can convert discy dwarfs 
to gas-poor spheroidal systems, but only if a heating source 
is imposed to keep the gas in the dwarf extended and hot 
at a temperature of ~ 2.5 x lO^K (L. Mayer, private com- 
munication). As we have shown, temperature floors of this 
magnitude inhibit the formation of thin discy dwarfs in the 
field. If a puffy dwarf of the kind we expect falls into the po- 
tential well of a larger galaxy with an extended hot gas halo, 
it should be quite susceptible to gas loss and morphological 
transformation. Trends of this kind do seem to be broadly 
in accord with some observations (e.g. Lisker et al. 2007, 
Geha et al. 2006). In another investigation, Mastropietro et 
al. (2005) used N-body simulations to study the transfor- 
mation of discy dwarfs to spheroidal dwarfs. They had some 
success, but were unable to reproduce the observed fraction 
of spheroidal dwarfs with negligible rotational support. Our 
results suggest that field dwarfs are likely to be born thicker 
than typically assumed and are therefore more susceptible 
to kinematic transformations. 

It has long been recognised that galaxy formation must 
become increasingly inefficient in dark matter haloes from 
the scale of big spirals to small dwarfs (White & Reese 
1978; Klypin et al. 1999; Moore et al. 1999; Strigari et al. 
2007). Our work suggests that it is difficult to avoid the 
suppression of galaxy formation efficiency in small haloes. 
Shallow potential wells naturally give rise to puffy galaxies 
with long star formation time scales. This situation makes 
them more susceptible to other feedback effects and exter- 
nal influences that may act to suppress star formation even 
further. Small galaxies are also systematically more metal 
poor than larger systems (e.g. Tremonti et al. 2004). The 
lower efficiency of star formation may also explain the ob- 
served mass-metallicity relation without the need for strong 
winds (Tassis et al. 2006). For example. Brooks et al. (2007) 
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recovered in their cosmological simulations the observed re- 
lation only by including an allowance for ISM heating (see 
also Ricotti & Gnedin 2005). 

A natural extension of this work would be a more 
detailed modeling of the coupling between various energy 
sources and the ISM. Ideally, this would include a star for- 
mation prescription based on molecular cooling and radia- 
tive transfer for the treatment of the feedback from stars 
to the ISM - allowing to overcome the simple proxy of a 
temperature floor. 
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